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We investigate the elementary excitations of quasi one- 
dimensional 5=| systems built up from zig-zag chains with 
general isotropic exchange constants, using exact (Lanczos) 
diagonalization for 24 spins and series expansions starting 
from the decoupled dimer limit. For the ideal one-dimensional 
zig-zag chain we discuss the systematic variation of the basic 
(magnon) triplet excitation with general exchange parameters 
and in particular the presence of practically flat dispersions in 
certain regions of phase space. We extend the dimer expan- 
sion in order to include the effects of 3D interactions on the 
spectra of weakly interacting zig-zag chains. In an application 
to KCUCI3 we show that this approach allows to determine 
the exchange interactions between individual pairs of spins 
from the spectra as determined in recent neutron scattering 
experiments. 



I. INTRODUCTION 

Spin systems consisting of chain- or ladderlike struc- 
tures as basic building blocks have recently attracted 
much attention. These systems are of interest on the 
one hand as one-dimensional (ID) model systems allow- 
ing to study quantum phase transitions related to the 
existence of a spin-gap and their dependence on the ex- 
change parametersElfl; on the other hand they describe an 
increasing number of real materials when an additional 
(small) exchange ooupling in the remaining two dimen- 
sions is introducecE. A material of particular recent ex- 
perimental interest is KCUCI3E 
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FIG. 1. Zig-zag chain as defined by ( |]J) 

We have performed an investigation of the dynamics 
of such systems with a twofold aim: (i) We discuss the 
spectrum ui(q x ) of the low- lying triplet excitations in the 
ideal ID system over a wide range of exchange parameters 
using both series expansions and exact diagonalization, 
in order to determine the range of applicability of the 
series expansion approach and to study the validity of 
using an effective interaction between dimers. We find 



and discuss in particular a regime in phase space with 
extremely small dispersion and a minimum of u>{q x ) at 
finite wavevector < q x < n. (ii) We extend the dimer 
expansion to include 3D couplings and apply this method 
in particular to a discussion of the dynamics of the quasi 
zigzag-ladder material KCuCbj in terms of microscopic 
exchange parameters. 

Of particular interest as an ID building block for this 
type of materials is the S = | zig-zag chain, as shown in 
Fig. [l] and defined by the following hamiltonian: 



H — JS\ : iS2,i + JiSi t iS2,i+i 

i=l 

+ J% f SijiSiji+l + S%,iS-i,i+\\ ■ (1) 

On the theoretical side, this generic ID model inter- 
polates between a number of seemingly different limiting 
cases: It is an alternative way to formulate the hamilto- 
nian for the generalized S = g spin ladder generalized 
to include one diagonal interaction or equivalently the 
5 = | chain with nearest neighbour (NN) alternating ex- 
change and next nearest neighbour exchange (we use the 
shorthand NNNA-chain in the following). It thus covers 
the well-known limiting models of the isotropic S = 5 
Heisenberg chain (HAF, J = Ji, J2 = 0), the standard 
antiferromagnetic S = h ladder (J\ = 0, J = J2 > 0), 
the weakly interacting dimer chain (Ji, J2 -C J) and 
the S = 1 antiferromagnetic (Haldane) chain (Ji —> 
—oo,J+ 2.J2 > 0). It can alternatively be considered 
as a two legged spin ladder with rung coupling </, leg 
coupling J2 and additional diagonal coupling J\. 

The theoretical interest in the dynamics of the NNNAp 
chain goes back to the work of Shastry and SutherlandEl, 
who identified the elementary excitations without alter- 
nation as free particles (spinons), which may become 
bound. A variational approach to the|-excitations of 
the NNNA-chain based on this concepfQ has recently 
been shown to cover qualitatively the transition from free 
spinons to the Haldane triplet. 

In a first approximation real materials are often consid- 
ered as examples of ID chains with a hamiltonian as given 
in Eq. (1); then they realize different points in the phase 
diagram spanned by the interaction constants J\ , J2 and 
illustrate the interest to describe systematically the vari- 
ation of static and dynamic properties with the parame- 
ters J\,J2- Static properties, such as susceptibility and 
specific heat, however, have turned out to be rather iru 
sensitive to the details of the raicxpscopic hamiltonianB 
and, in the case of (VO)2P2 07lTt!J have even not been 
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able to reveal the basic interactions as two- instead of 
one-dimensional. 

Thus for a description of real materials a systematic 
microscopic treatment of the dynamics is of particu- 
lar importance. In section II we present a systematic 
overview of the dynamics of the ID system, discussing 
both general properties as well as comparing results from 
exact diagonalization to results from series expansions. 
In section III we use the series expansion approach to cal- 
culate the low-lying excitations in a 3D material with the 
structure of KCuCLj in terms of the microscopic hamilto- 
nian. This will allow us to go beyond the determinates 
of effective dimer exchange parameters in recent worker 
and to determine the microscopic exchange parameters 
A summary will be given in the concluding section IV. 



II. ELEMENTARY EXCITATIONS OF THE ID 
ZIG-ZAG CHAIN 

We start with a short summary of the symmetries of 
the ID system (the chain direction is denoted as x— axis): 
Translational symmetry is described by the wavevector q x 
defined in a Brillouin zone — — < q x < + -; the unit cell 
of length a contains two spins, or equivalently one dimer 
(one singlet in the limit J± = J 2 = 0). We thus use 
the conventional notation for ladders: a is the distance 
between rungs, whereas the distance between spins in the 
NNNA-chain picture is a/2. We thus expect two basic 
excitations per unit cell. We will use units a = 1 in the 
following. Excitation frequencies at wavevectors q x and 
— q x are equal owing to reflection symmetry along the 
chain. 

For special points in the phase diagram additional sym- 
metries exist: 

Without alternation (Ji = J) it is natural to use a 
unit cell of length a = containing only one spin. Our 
Brillouin zone is half of this Brillouin zone of the uniform 
chain and the excitations of the conventional spin chain 
will appear folded back to our smaller Brillouin zone. 

For the ladder symmetry ( J\ — 0) there exists a quan- 
tum number parity, P, resulting from the interchange 
of the two legs and we can classify states as positive or 
negative under this reflection. Each dimer in the singlet 
(triplet) state contributes a factor of — 1(+1) to this par- 
ity. An alternative notation introduces the component 
q± with values (corresponding to P = +1) and 7r (cor- 
responding to P = — 1). 

The ground state is a singlet in the whole phase plane 
and the lowest excited state is generally a triplet. The 
ground states for the ladder symmetry (Ji = 0) have 
parity P = +1 for L even. 

We have studied the dispersion Lo(q x ) of the basic 
triplet excitation for a typical variety of paths in the 
J\ — J2 parameter space by two methods: (i) By exact 
numerical diagonalization, using the Lanczos algorithm, 
we have calculated cu(q x ) for the lowest excited states 



(between 2 and 4 states) for 24 spins, i.e. for 7 different 
values of wavevector q x . (ii) We have performed series 
expansions around the dimer point, J\ — Ji — 0, up to 
third order analytically and up to 10th.|-Qrder after im- 
plementation of the cluster algorithmtHlli^l on an Alpha 
work station. Thus we have obtained the ground state 
energy Eq and an effective Hamiltonian which can be di- 
agonalized by a Fourier transformation. Finally we get 
the dispersion relation for the lowest excited state ex- 
pressed as series J2 n a ™ cos(nq x ). 

In the follwoing we present a number of results for the 
ID zig-zag chain which prepare the stage for the first 
application of the method to a nontrivial 3D system in 
section III and also add some new aspects to the large 
number of previous studies on the ID system defined by 
eq. (1) in recent years. To give a short review of existing 
work we mention first that the dimer series expansion 
approach startecLwhen the work of Brooks Harriaiil was 
revived by Uhrigij in the context of CuGe03. The ex- 
pansion for the triplet disperskp was extended to high 
orders recently by Oitmaa et aO for ladders to 8 th order, 
by Barnes et altZl for the Heisenberg alternating chain to 
9 th order and by Singh et alia for the disorder line to 23 th 
order (using a special symmetry on this line) . The model 
of Eq. (1) was also treated Jby alternative methods as 
random phase approximationtH Bruckner theory for the 
equivalent dilute Bose gaaHlrcil, exact diagonalizationEL 
and continued fraction expansion based on ED resultsEEl 
and DMRGcj. From these studies a rather complete pic- 
ture of the low-energy dynamics of the ID zigzag chain 
has emerged. In this section we supplement this picture 
by the following two remarks. 

We start from the neighbourhood of the dimer point 
where an expansion in J\, Ji to low orders is sufficient. 
Up to third order the following result for the dispersion is 
obtained (as given in ref. |l^, frequency and exchange con- 
stants are measured in units of the intradimer exchange 
J from now on): 
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In first order perturbation theory in J\ and Ji the 
spectrum is dispersionless on the Shastry-Sutherland line 
Ji = 2J2 (also known as 'disorder line'). Eq. (2) also 
shows the general feature that the location of the mini- 
mum of the dispersion curve shifts from q x = to q x = tt 
somewhere close to crossing this line. 

In recent approximpt»,theoretical treatments of inter- 
acting dimer systemseJO the exchange interactions be- 



tween individual were reduced to an effective interaction 
between dimers. For the present ID zigzag chain the re- 
sult of this approximation is that only the combination 
J2 — i Ji enters into the dispersion. In more detail, the 
result is 



w{q x ) = Jl + 26u>W{q x ) 



(3) 



where 1 + 5a/ 1 ) ((fc) is the dispersion in lowest order, 
i.e. the result of simple propagation of an excited dimcr 
triplet without considering its coupling to higher energy 
modes. Evidently this is true in lowest order of the ex- 
pansion; it is seen, however, already from Eq. (2) that ad- 
ditional terms which depend on the individual exchange 
interactions enter in higher order. Comparing the higher 
order coefficients in the series expansion we find that the 
effective dimer approximation of Eq. (3) amounts to keep- 
ing only the leading (i.e. lowest) powers in J\, J2 for each 
coefficient a n of cos(nq x ). 




FIG. 3. Variation of the dispersion curves when the in- 
commensurate region is crossed. Inset shows with high resolu- 
tion an example of a dispersion curve with minimum at finite 
wavevector as obtained from exact (Lanczos) diagonalization 
as well as from the dimer expansion method. 
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FIG. 2. Comparison of excitation spectra for fixed effec- 



tive dimer interaction J2 — i J\ — 



FIG. 4. As Fig. but for different couplings 



For a quantitative check of the effective dimer approxi- 
mation we refer first to ref . [l8| where the dispersion on the 
disorder line, i.e. for fixed J 2 — \J\ = is shown. The 
effective dimer approximation is reasonable for a large 
part of the line but deteriorates rapidly when the non- 
alternating limit Ji = 1, i.e. the Majumdar- Ghosh point 
is approached. For the more generic value J2 — 5 J\ = | 
we show the dispersion in Fig. ||; it is seen that the ef- 
fective dimer approximation is of very limited value in 
this case. In addition these results, comparing the series 
expansion results to spectra from exact diagonalization, 
we find the following limits of validity for the dimer series 
expansion method: Whereas the 10 th order spectra coin- 
cide with exact diagonalization results within a few per- 
cent for |Ji| < 1 and J2 < 0.5 the accuracy deteriorates 
rapidly when the symmetric ladder ladder configuration 
(Ji = 0, Ji = 1) is approached. 



Whereas the dispersion is flat on the disorder line to 
first order in J\ , J 2 , but develops maxima and minima 
in higher orders with an apparent jump of the minimum 
energy form q x = n to q x = tt close to this line, it was 
noted already beforeEm from exact diagonalization results 
that a regime with extremely flat dispersion and a min- 
imum at finite wavevector < q x < n might exist in a 
narrow regime in J\ — J 2 parameter space. We have in- 
vestigated this point again using the series expansions to 
10 th order which provide us with a continuous wavevector 
dependence and have confirmed the earlier speculation: 
In Fig. |^, |] we show two examples for spectra on the 
lines Ji = 1 — 2J 2 resp. J\ = Ji)'- The points in 

Fig. U were discussed before in ref.H; Fig. || demonstrates 
that the shallow minimum of the dispersion curve at a 
wavevector in the middle of the Brillouin zone is repro- 
duced in the series expansions. As further example we 
show in Fig. H the flat dispersion for somewhat smaller 



values of J\. In Fig. |^ we show the development of the 
wavevector g m i n for the minimum value of u>(q) on the 
line Ji = |(1 — J 2) in the J\ — J2 parameter space (cor- 
responding to the case of Fig. ||). The gradual transition 
of <7min between g m i n = and g m i„ = tt, possibly with 
infinite slope at the end points, is clearly seen. The com- 
parison between results in 8 th and 10 th order illustrate 
the convergence of the expansion. At present, however, 
the series expansions do not give a hint to a possible 
fundamental reason for nearly flat dispersion: We have 
examined the expansion coefficients a n of cos(nq x ) up to 
10 th order, but we do not find any indication of a con- 
vergence to zero for n ^ 0. 



• 10th order 
X 8th order 



0.2350 0.2355 
J, 



FIG. 5. Wavevector g m i n for the minimum value of u>(q) 
on the line Ji = |(1 — J 2) 
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FIG. 6. KCUCI3— structure projected on the xz— plane. 
Thick full lines denote the basic dimers, the height above the 
reference plane y = nb is given for the two spins of one dimer 
in each zig-zag chain. Thin lines show interdimer interac- 
tions: Full thin lines denote interactions of spins in dimers 
with identical values of n, dashed and dotted thin lines de- 
note interactions of spins in dimers with different values of n. 
As indicated by the arrows, the dashed (dotted) lines start at 
n and end at n — 1 (n + 1). Only one interaction of each type 
is shown. 



III. INTERACTING ZIG-ZAG CHAINS 

In real materials consisting ofrWjeakly interacting chains 
such as KCuCyaEl and CuGeOsES it is clear from inelastic 
neutron scattering experiments that there is considerable 
dispersion for wavevectors perpendicular to the double 
chain direction. For weakly interacting chains, the series 
expansion approach is to be considered as the only reli- 
able systematic approach which (by comparing results in 
subsequent orders) allows a consistency check. Whereas 
the series expansions have been extended to cover sys- 
tems coupled in 2D and interesting results have been ob- 
tained for the spin Peierls materiaLCuGeOsEj, for the 1/5 
depleted scumre material CaV^gEIl as well as for general 
parameters^, we present in the following the first results 
using the expansion for weakly interacting dimers for a 
3D coupled system with particular application to inves- 
tigate the magnon dispersion in the material KCUCI3. 
The zig-zag chain system KCUCI3 actually appears to be 
closer to the dimer point than the systems treated so far 
and therefore is supposed to be a better candidate for 
this expansion. 



The structure of KCUCI3 is shown schematically in 
Fig. ^| in a projection to what is conventionally called the 
xz— plane: The fundamental dimers which are shown as 
solid lines form zig-zag chains in the x— direction, neigh- 
bouring zig-zag chains are shifted with respect to each 
other in y— direction by half a lattice constant |; this 
shift as well as the tilting of the internal dimer direction 
are indicated in Fig. |^ by giving the y— coordinate for 
each line of spins in a;— direction (n is an integer which 
numbers the different planes). The elementary cell con- 
sists of two dimers, dimer D\ at the origin Ri = and 
dimer D2 at position R2 — (0, 5, 5) and the two spins 
forming each dimer, 1 and 2, are at positions Ri + di for 
spin 1 and Ri — di for spin 2, for dcfinitcness we take 
<k,z > 0. 

For the exchange interactions between spins we will use 
the following notation: The main intradimer exchange 
is denoted as J. The exchange interaction per bond 
between spins in dimers separated by a lattice vector 
R = lae x + mbe y + nce z is denoted as J(i rnn ) for the 
exchange between equivalent spins (pairs (11) or (22)) of 
the corresponding dimers and as J'n mn \ for the exchange 
between nonequivalent spins (pairs (12) or (21)). The 
following exchange interactions will be considered: J(ioo) 
for pairs (11) and (22), J/ 10 ol ^ or ^ ne P a ^ r (-^) ano - "^(201) 



for the pair (12), and for the two cases p = and p = 1: 
<^( p ii) fo r the pair (11) starting from dimer D\ and for 
the pair (22) starting form dimer D2 , J( P -^ I) = ^(p||) 
for the pair (22) starting from dimer Z?i and for the pair 
(11) starting from dimer D2, J', 



(p- 



J', 

(p 



n> for pairs 

22/ 



(12) starting from either dimer D\ or dimer Z?2- 

Alternatively we can look at the structure projecting 
on the plane spanned by the directions e y and e x + ^e z . In 
a schematic picture which shows the topology of the ex- 
change interactions only, we obtain Fig. 0; the 3D struc- 
ture of KCUCI3 results when identical planes are stacked 
and the fundamental dimers are connected by zig-zag in- 
teractions. Evidently the planar structure of Fig. [?] can 
be reduced to a number of limiting cases including cou- 
pled alternating chains (J^ 1± ii^ = 0) or coupled ladders 

(^(201) = ^u, = 0). 

If the dispersion is considered only to first ordero 
or if higher .-.orders are included in an RPA-like 
approximation^, J and J' for a given (Imn) enter only in 
the combination J cff = ^(pJ — p'J'), where p,p' = 1,2 
is the number of exchange paths between the interacting 
dimers; J off is denoted as effective dimer interaction. As 
already seen in the ID case of section 2, a correct treat- 
ment beyond first order involves J and J' independently. 
According to previous workcm the main exchange interac- 
tions in addition to the basic intradimer exchange are be- 
tween dimers separated by (Imn) = (100), (201), (l±^i). 
Except for (201) these dimer-dimer interactions involve 
both J and J' and it is our aim in the following to dis- 
cuss the validity of the effective dimer approximation for 
KCUCI3 and to investigate to what extent the exchange 
interactions between individual spins can be determined 
from the present status of experimental results. 

A calculation up to second order in the ratios J(i m n) I J 
leads to the following expression for the frequency of the 
basic triplet (frequency and coupling constants are mea- 
sured in units of the basic dimer exchange constant J 
and wavevectors qi are given in units with the crystallo- 
graphic lattice constants a, b, c set equal to unity): 

u(q) = l + SuJ W (q) - ^Sj 1)2 (q) 

+ J (100) (j(100) _ ^(100)) - ^^'(100) COS ^ 



+ J (°H) J (oH) + J (iH) J (iH) 



2 i J (Wir J '(°H) 



q v q z 

COS — COS — 

2 2 



+ 2{ J (^r J (^)) cos Y cos - 



l 2 

^'(201) COS ( 2 9z 



(4) 



Here 



ScjW(q)= i (2J (100) - J ( ' 100) ) 



cosq x 



J, 



J (°U)) 



cos -f cos— 
2 2 



(our j ( 

J H!r J (iH)J cos T cos ^^ 

-J ( ' 201) cos(2 fc + q z ) 



(5) 



is the dispersion to first order in the exchange constants, 
i. e. the effect of simple propagation of an excited dimer 
triplet. 

In Eq. |5| we have used for simplicity an extended zone 
scheme: Since there are two dimers in the proper crystal- 
lographic unit cell, there are two branches of triplet exci- 
tations for each wave vector in the crystallographic Bril- 
louin zone (— tt < q y b, q z c < 7r). Owing to the symmetry 
of the exchange interactions in the hamiltonian, these two 
branches join to the unique smooth expression given in 
Eq. [|, when the doubled zone 2n < q z c + q y b, q z c — q y b < 
2ir is used. Each triplet excitation, however, is present 
at all equivalent wavevectors of the crystallographic re- 
ciprocal lattice. For a discussion of dipole transition am- 
plitudes the noninteracting triplet approximation can be 
used as sufficient guide and gives the follwing resultsO: 
For momentum tranfer perpendicular (parallel) to e y the 
branch of Eq. || in the first (second) crystallographic Bril- 
louin zone gives the only nonvanishing contribution. (A 
change form the first to second crystallographic Brillouin 
zone then corresponds to a change in either q y or q z by 
2ir, i.e. to a change of sign in the second and third line 
ofEq. |). 




1 (20 1) 



J (1 1/2 1/2) J (1 1/2 1/2) 

FIG. 7. Alternative schematic view on the structure of 
KCuCl 3 . 

We note that the wavevector q z should be distinguished 
from the quantity q± which is often used to denote the 
two values of the quantum number parity discussed in 
section II as q± — 0, ir. This wavevector is measured in 
units of c _1 , where c is the rung length, which differs 



from the reduced lattice constant |. In experiments so 
far only the basic triplet with q± = it has been observed 
and the interchange of minima with variation of q z is an 
effect of the crystallographic lattice geometry and not of 
the ID ladder geometry. Excitations with q± — are 
excitations with two dimer quanta and have a minimum 
energy of twice the gap energy A. 

In order to obtain results which are quantitatively re- 
liable we have performed series expansions for the 3D 
coupled system to 4 th order following the lines described 
in section II. Because of the complex lattice we did not 
characterize the clusters which results in a large number 
of clusters: 5532 clusters in 4th order. The convergence 
of these series expansion results is excellent for the small 
values of the expansion parameters J(i mn ) /J < 0.4 in the 
application to KCUCI3. 
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FIG. 8. Dispersion curves for KCUCI3 in various directions 
in q— space showing the variation with J'. 1 ± . at fixed effective 

dimer interaction J,, 1 i\ — J,' 1 1, 

1 22' ^iiJ 
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In Figs. H and § we show results for dispersions along 
typical lines in if— space together with the data points 
from ref. |[ Dispersions are plotted only for the branch 
which has nonvanishing dipole transition amplitude in 



lowest order leading to + or — sign in Eq. || as dis- 
cussed above. We assume that the effective dimer 
change takes the values determined in previous workl 
and have therefore fixed the following combinations of 
exchange parameters: 
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= 0.188, 
-0.110, 

0.160 (6) 




J' (100) =0.300 
- J'aoo,=Q-500 



FIG. 9. Dispersion curves for KCUCI3 in various directions 
in q— space showing the variation with J( 100 ) at fixed effective 
dimer interaction J100 — 5^/(100) ~ —0.055 

The interactions J( ,i,i): J[q i i) are found to be negli- 
gibly small. In order to demonstrate the relevance of the 
individual spin exchange parameters (as opposed to the 
effective description) we discuss the separations (100) and 
(l±^i) independently: In Fig. || we present the variation 
of the dispersions w(q) in selected directions in q space 
for different values of the coupling J(\ii) = 0,0.2,0.4 
with the remaining parameters fixed as given above. It 
is seen that the distribution of the effective dimer in- 
teraction between parallel and diagonal terms essentially 



shifts the dispersion curve by constant amounts. A com- 
parison to the corresponding neutron scattering results 
leads to the conclusion that J(iii) = 0.200 is the most 
likely value. The analogous results for different values of 
the diagonal (zig-zag) coupling J/ 100 ) = 0.1,0.3,0.5 are 
shown in Figs. ^|. Here the frequency ui(q x = 0) (which is 
the energy gap for q z = tt and the dip energy for q z = 0) 
depends only on the effective interaction, whereas the 
frequency u)(q x = tt) (which is minimum for q z — and a 
dip energy for q z = 2tt) allows to determine the exchange 
between individual spins. 

Comparison of our series expansion results to the neu- 
tron scattering data then leads to the following values for 
the exchange constants between individual spins, not yet 
determined by the values for the effective dimer interac- 
tions published so far: 



It is shown that the effective dimer approximation, when 
treated in the random phase approximation, sums up the 
leading powers of the dimer expansion. 
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(100) 



-0.005 



J ( ' 100) « 0.100. 

These values imply that the ladder system in KCUCI3 
is much closer to an alternating spin chain than was be- 
lieved so far and the leg interaction tends to be ferromag- 
netic if it is nonzero at all. The results for the interchain 
interactions in (I55) direction are much less conclusive. 
The most likely values are 

J(iH) w °- 200 

X'ih« 0.040. 

However, the error is large and data actually may be 
compatible also with smaller values for 



(Hi)- 



IV. CONCLUSIONS 

We have investigated the dispersion curve for the 
low energy triplet excitations of one-dimensional and of 
weakly coupled zig-zag chains starting from the limit of 
nonintcracting dimers and performing an expansion in 
the interdimer interactions. The series up to 10 th order in 
the ID case and up to 4 th order in the 3D case were eval- 
uated explicitly after implementation on a work station. 
In the ID case the dispersion curves agree with those 
obtained from exact diagonalization using the Lanczos 
algorithm for a large regime around the dimer point; in 
this regime the method provides a reliable approach to 
calculate the dispersion in its continuous dependence on 
the wavevector. In a narrow regime close to the Shastry- 
Suthcrland line we find the minimum of the dispersion 
curve at an intermediate wavevector fc m i n , < k m i n < tt. 
As an application of the 3D case we present the applica- 
tion of the method to the KCuCbj — structure. Fitting to 
the dispersion as measured in inelastic neutron scatter- 
ing experimentsEl we determine the exchange interchange 
interactions between individual spins in addition to the 
effective interaction between dimers determined before. 
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